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ABSTRACT 


The standard tropospheric calibration model implemented in the operational 
Orbit Determination Program is the seasonal model developed by C. C. Chao in 
the early 1970s. The seasonal model has seen only slight modification since its 
release, particularly in the format and content of the zenith delay calibrations. 
Chao’s most recent standard mapping tables, which are used to project the zenith 
delay calibrations along the station-to-spacecraft line of sight, have not been 
modified since they were first published in late 1972. This report focuses 
principally on proposed upgrades to the zenith delay mapping process, although 
modeling improvements to the zenith delay calibration process are also discussed. 
A number of candidate approximation models for the tropospheric mapping are 
evaluated, including the semi-analytic mapping function of Lanyi, and the semi- 
empirical mapping functions of Davis et al. (“CfA-2.2”), of Ifadis (global solution 
model), of Herring (“MTT”), and of Niell (“NMF”). All of the candidate 
mapping functions are superior to the Chao standard mapping tables and 
approximation formulas when evaluated against the current Deep Space Network 
Mark III intercontinental very long baseline interferometry database. 
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SECTION 1 


INTRODUCTION 


An integral part of NASA’s Deep Space Network (DSN)-based tracking and navigation 
system is the collection of ground-based observations, which are used to continuously correct the 
prediction of a spacecraft’s position in time. These observations, made by a global network of 
DSN radio antennas, consist of radio metric tracking data such as: 1) Doppler, a measure of 
“frequency change” between the spacecraft and ground station; 2) range, a measure of line-of- 
sight “distance” between the spacecraft and ground station; and 3) Very Long Baseline 
Interferometry (VLBI), an interferometric measure of the spacecraft’s “angular position” on the 
plane of the sky. Radio signals that pass through the Earth’s atmosphere are sensitive to 
electrically neutral particles which cause a time delay in the signal, depending on the dielectric 
constant of the media and path length. Additionally, if the signal does not travel parallel to the 
density gradient of the atmosphere, it will be subject to atmospheric refraction, causing its path to 
deviate from a straight line. Geometrical bending of the signal path is greatly intensified in the 
lower atmosphere, most notably the troposphere, due to the curvature of the air layer near the 
Earth’s surface (see fig. 1). This so-called “bending effect” is significant at low station-to- 
spacecraft elevation angles (< 10°). Consequently, improved modeling techniques for the excess 
propagation delay due to the neutral atmosphere must be incorporated into the operational 
navigation software in order to take full advantage of the high-fidelity zenith delay calibrations 
that are provided by the DSN’s Tracking System Analytic Calibration (TSAC) activity, and 
further, to reduce the effect of elevation-dependent “mapping” errors induced by projecting the 
calibrated zenith delays along the station-to-spacecraft line of sight. In this report, the refraction 
and zenith delay effects of the entire neutral atmosphere are collectively referred to as the 
“tropospheric effect.” 

The tropospheric calibration model currently implemented in the operational Orbit 
Determination Program (ODP) 1 of the Jet Propulsion Laboratory (JPL) is based in large part on 
the seasonal model developed by C. C. Chao of the Mariner Mars 1971 (MM 71) TSAC activity 
[Chao, 1974], Additionally, in past and current planetary missions, two-way radio metric 
tracking data acquired from Deep Space Stations (DSSs), have frequently been either deleted (cut 
off) below 10° to 15° elevation angle or artificially “deweighted” by an elevation-dependent 
empirical weighting function. The weighting function was developed as an attempt to alleviate 
any modeling errors introduced by the media calibration process. Unfortunately, these ad hoc 
procedures have pre-empted the need to upgrade the media calibration models in the ODP. The 
argument for eliminating low-elevation data clearly breaks down when differential data types are 
incorporated into the orbit determination process; relatively low elevation angles (~6°) must be 
used in the differential data acquisition process in order to obtain sufficient tracking coverage 
between co-observing stations on DSN intercontinental baselines. Differential data types such as 


I The ODP is a large institutional software system used for research and navigation support of flight operations 
[Panagiotacopulos et al 1974]. 




Figure 1. Geometrical bending of the signal path due to tropospheric refraction. 

Note: E 0 represents the observed, or “refracted,” elevation angle, while E denotes the 
“unrefracted” (straight-line) elevation angle. (The size of the tropospheric air layer and ray 
bending are exaggerated for illustration purposes.) 

wideband and narrowband VLBI and differenced Doppler and ranging are often used to provide 
a direct angular measure of the spacecraft on the plane of the sky, as they are typically less 
sensitive to mismodelings of spacecraft nongravitational forces than conventional two-way radio 
metric data types. 

In Section 2 of this report, the ODP’s currently implemented tropospheric calibration model 
is extensively reviewed. Because this is intended to be an archival document, where appropriate, 
all relevant mathematical formulations are provided. Section 3 discusses a number of proposed 
modeling improvements in both the zenith delay calibration and zenith delay mapping processes. 
Section 4 lays the mathematical framework for a number of candidate tropospheric mapping 
approximation models, commonly referred to as “tropospheric mapping functions,” which has 
been a very active area of research during the past 2 decades (see e.g., Gallini [1994]). Section 5 
presents a statistical evaluation of the candidate tropospheric mapping functions based on an 
extensive set of VLBI measurements acquired at the DSN over a 5 year period extending from 
1988 to 1993 using the Mark III data acquisition systems. Section 6 discusses general selection 
criteria and implementation issues of the candidate mapping functions. Section 7 concludes the 
report with a general summary statement and specific recommendations for modernizing the 
tropospheric calibration model in the ODP. 
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SECTION 2 


REVIEW OF THE CURRENT MODEL 


An excellent discussion of the transition from early ODP tropospheric calibration models to 
the currently implemented seasonal model is given by Mottinger [1984]. There have been a few 
changes to the standard model since the time of Mottinger’ s report, which are documented here 
along with a review of the current implementation process. In the ODP, the seasonal 
tropospheric calibration model comprises two fundamental parts: 1) one-way range zenith delay 
calibrations, or “adjusts,” supplied by TSAC in the form of coefficients of a normalized 
polynomial or Fourier series in time, and 2) standard mapping tables, which are used to project 
the zenith range correction to the direction of the source, as viewed by the observer. The source 
can be any spacecraft or extragalactic radio source (EGRS), i.e., quasar, and the observer can be 
any DSS. A compact formula is then used to compute the actual range tropospheric correction 
along the line of sight and convert the result to appropriate units, consistent with the data type 
being used ( e.g ., Doppler, range, VLBI). This entire process is handled in the program link 
EDIT independently from the light-time determination and observable formulation, which is 
performed in the program link REGRES (see Fig. 2) [Ekelund, 1993]. The standard 
implementation makes use of the common practice of separately modeling the hydrostatic, or 
“dry,” component (which comprises atmospheric constituents in hydrostatic equilibrium- 
principally, N 2 , 0 2 , Ar, and C0 2 ) from the highly variable water-vapor, or “wet,” component. 
Fortunately, the major contributor to the total delay is the dry component (-90% to -95%), and 
its variability at any DSS site is small and seasonal in nature. Typical observed values at the 
Deep Space Communications Complexes (DSCCs) of DSCC 10 (Goldstone, California), DSCC 
40 (Canberra, Australia), and DSCC 60 (Madrid, Spain) for the zenith dry and wet delays are 
-200 cm to -216 cm and -1 cm to -25 cm, respectively [Runge, 1993]. 

2.1 Zenith Range Adjusts 

The TSAC activity employs a seasonal model to produce calibrations for the zenith dry and 
wet tropospheric delays and delivers them to flight project navigation teams as coefficients of a 
normalized polynomial or Fourier series in time. The original seasonal model was derived by 
Chao in the early 1970s, during which time he studied the influence of seasonal fluctuations in 
surface meteorological measurements on refractivity, based on radiosonde balloon measurements 
acquired over a 2 year period (1967 and 1968) from six globally distributed weather stations 
[Chao, 1971]. Chao concluded that the refractivity profiles repeated satisfactorily over the 2 year 
period and that monthly averages of surface measurements of four meteorological parameters, 
i.e., surface barometric pressure, P 0 , linearly extrapolated temperature, To, temperature lapse rate 
(altitude temperature gradient), w, and surface relative humidity, (RH)o, could be used to predict 
the zenith dry and wet one-way range delays within a reasonable degree of accuracy, provided 
that corrections were made to Po and To to account for any station altitude offset between the 
radiosonde station and the nearest DSS. Zenith delay calibrations were determined from polyno- 
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Figure 2. ODP software set functional diagram for observable generation, data editing, and 
accumulation of the information matrix. 

Note: The program links that contain the media models reside in program links REGRES and EDIT. 
Also shown are the ODP input/output file interfaces. The file extension .NIO represent direct access, 
binary NAVigation Input-Output (NAVIO) files, while the file extension .TXT represents text-based 
input files, generally assumed to be Fortran NAMELIST files. 


mial fits to zenith delay values computed from monthly averages of surface meteorological 
measurements using the “Berman 70” model [ Berman , 1970], which was derived under the 
assumption of static equilibrium, perfect gas law for the troposphere, and constant gravitational 
acceleration: 
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A = 17.1485 

B = 4684.1 

C = 38.45 

g = gravitational acceleration, 980.6 cm/s 2 
R d = perfect gas constant of dry air, 0.287 J/g*K 

g/R d = 34. 1 K/km 

C, = 77.6 

C 2 = 29341.0 

In the seasonal model, monthly means of ( RH) 0 were averaged over each month to overcome the 
inadequacy of Eq. (2) to possible fluctuations in relative humidity with altitude. This model was 
adopted as the standard zenith delay tropospheric calibration method for MM'71 [Chao, 1974], 

Although the format and content of the TSAC tropospheric calibration polynomials has 
evolved over the years, they are still delivered as “card images” which conform to the Command 
Statement Processor (CSP) English-like command language used by the ODP to adjust, delete, or 
weight data points. 2 In general, a CSP command consists of a verb and optional scope limiters 
and computation specifiers, terminated by a period [DPTRAJ-ODP User’s Reference Manual, 
Vol. 1, Rev. F, 1994], A command may, and frequently does, extend over more than one card 
image. Elements of CSP commands used in the TSAC media calibrations are summarized in 
Table I. 

The ADJUST command for media calibrations takes the general form 

ADJUST (data-type) BY SERIES-TYPE (coefficients) 

MODEL (model-name) FROM (A) TO (B) DSN (complex) 

Additional explanatory information, not to be used in the actual adjustment process, may appear 
on a command provided that it is preceded by a number sign (#). The ADJUST (data-type) 
command verb indicates that the calibrations were derived especially for a specific class of data 
types. The allowable options for the data-type argument are ALL (corrections to be applied to all 
data types), DOPRNG (Doppler and range), or VLBI (exclusively). 

Media calibrations may be specified by a normalized power series (NRMPOW), a Fourier 
series (TRIG), or in special cases, by a single constant (CONST). In the past, the coefficients or 
the constant could be represented as single- or double-precision floating-point numbers with the 
computation specifier preceded by the letter “D” for the double-precision format ( e.g ., 
DNRMPOW). This level of precision has since been determined to be unnecessary, hence, all 

•j 

A detailed description of the format and content of TSAC media calibrations is given in DSN Tracking System Interfaces , Media 
Calibration Interface (TRK-2-23) [1985]. Excerpts from that document are given here for completeness and include recent 
undocumented modifications and enhancements. 
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coefficients or constants are provided as single-precision floating point numbers accurate to four 
decimal places [Royden, 1994], 


Table 1. Command Statement Processor (CSP) elements used in TSAC media calibrations. 


T yp e 

Verbs 


of Element 


Elements Used 


ADJUST, DELETE 


Scope Limiters 


Time Span Limiters: 

FROM BEFORE 

TO AFTER 

AT 


Network (Complex or Station): 

DSN 

CIO (Goldstone) 11,12, 19 

C40 (Australia) 41, 42, . . . , 49 

C60 (Spain) 61,62, ..., 69 


Band: S, X 


Source: SCID, QUASAR 


Computation Specifier 


Series Type: 
BY NRMPOW 
BY CONST 
BY TRIG 


Model: 

DRY NUPART (dry troposphere) 

WET NUPART (tropospheric water vapor) 
CHPART (ionosphere) 

DRVID (ionosphere + solar plasma) 3 

Data Type: 

ALL 

DOPRNG 

VLBI 


“obsolete 


For tropospheric calibrations, the MODEL (model-name) computation specifier takes the 
form of either MODEL (DRY NUPART) or MODEL (WET NUPART), depending on whether 
dry or wet zenith tropospheric delay calibrations are to be calculated from the polynomial 
coefficients supplied. Series representation of the zenith range adjusts are determined as follows: 
Let T be the time associated with the current data point. For the BY SERIES (C„, C,, ..., C„) 
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FROM (A) TO (B) part of the ADJUST command, the zenith range adjust, A pz, is computed 
from one of the following series types: 

• Nonnalized power series (BY NRMPOW) 


n 


&Pz = X C k X k 
k= 0 


(3) 


This is referred to as a “normalized” power series because the time argument of the 
series, T, is replaced by a normalized, dimensionless argument X, which is defined to 
be -1 when T - A and +1 when T = B. The normalization is performed through the 
change of variable X = 2((T - A)/(B - A)) - 1. Currently, TSAC uses a lOth-degree 
polynomial to perform the fit. Thus, for the standard model, n = 10. 

• Fourier series (BY TRIG) 

Renaming the coefficients as BY TRIG (P,C,D u E x< D 2 ,E 2 ,...,D n ,E„)\ 


Apz = C + ^ [D/c cos (2k k X) + 

*=' (4) 

Eic sin (2k k X)\ 

where P is the period of the fundamental mode (in seconds) and X - (T- A) / P. For 
the standard model, n = 4. 

In some cases, an additional adjustment may be given in the form of a simple constant 
such as a correction for a station altitude offset: 

• Constant (BY CONST) 


A pz = Co (5) 

The format of the time sequence is taken to be YY/MM/DD,HH,mm,ss.ff, where YY is the 
last two digits of the year, MM is the number of the month, DD is the day of the month, HH is 
the hour, mm is the minute and ss.ff the seconds and fractional part of a second, respectively. Up 
to three of the right-most time fields may be omitted. When omitted, their default values are 
taken to be zeros. 

The dry and wet adjust cards are repeated pairwise as necessary to cover all DSCCs from 
which radio metric tracking data are expected. These are denoted by the scope limiters DSN 
(CIO), DSN (C40), and DSN (C60). Two examples of tropospheric adjust card image streams 
are shown in Figs. 3a and 3b and are representative of the “official” seasonal model. Illustrated 
are the lOth-degree normalized power series representation (NRMPOW) and 4th-degree Fourier 
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# DSN TROPOSPHERE CALIBRATIONS FROM 89/1/1 TO 91/1/1 

# FOR MAGELLAN ODP 

# 

ADJUST (ALL) BY NRMPOW ( 0.0539, -0.1545, 0.6120, 0.8394, -1.7571, 

-1.1924, 1.7748, 0.6266, -0.7482, -0.1117, 0.1137) MODEL (WET NUPART) 

FROM (89/1/1,00:00:00.01) TO (91/1/1) DSN (CIO). 


ADJUST (ALL) BY NRMPOW ( 2.0601, 0.0000, -0.1453, -0.0003, 0.4194, 

-0.0010, -0.4272, 0.0029, 0.1822, -0.0014, -0.0281) MODEL (DRY NUPART) 

FROM (89/1/1,00:00:00.01) TO (91/1/1) DSN (CIO). 


ADJUST (ALL) BY NRMPOW ( 0.1393, 0.0148, -0.4576, -0.0854, 1.3414, 

0.1292, -1.3936, -0.0714, 0.6081, 0.0131, -0.0961) MODEL (WET NUPART) 

FROM (89/1/1,00:00:00.01) TO (91/1/1) DSN (C40). 


ADJUST (ALL) BY NRMPOW ( 2.1541, -0.0013, 0.0777, 0.0063, -0.2608, 

-0.0077, 0.3180, 0.0032, -0.1643, -0.0004, 0.0305) MODEL (DRY NUPART) 

FROM (89/1/1,00:00:00.01) TO (91/1/1) DSN (C40). 


ADJUST (ALL) BY NRMPOW ( 0.0986, -0.1422, 0.4960, 0.7966, -1.4194, 

-1.1690, 1.4176, 0.6323, -0.5845, -0.1150, 0.0858) MODEL (WET NUPART) 

FROM (89/1/1,00:00:00.01) TO (91/1/1) DSN (C60). 


ADJUST (ALL) BY NRMPOW ( 2.1128, -0.0014, -0.0614, 0.0021, 0.1593 

0.0076, -0.1330, -0.0108, 0.0392, 0.0033, -0.0026) MODEL (DRY NUPART) 

FROM (89/1/1,00:00:00.01) TO (91/1/1) DSN (C60). 

# 


Figure 3a. Format of the TSAC-supplied zenith delay tropospheric calibrations for the 
“official” seasonal model given as a lOth-degree normalized polynomial. 

Note: Sample time span from the Magellan mission. 


series representation (TRIG), respectively. Either one of these representations can be used for 
operational orbit determination; however, the latter has been shown to give a better 
characterization of the fitted data [. Royden , 1994], A further advantage of the TRIG 
representation is that the user does not need to change the time sequence for which the 
tropospheric corrections are to be applied, as is currently required in the NRMPOW case (which 
only covers a 2 year span). 

In the ODP, all ADJUST commands compute a correction which is subtracted from the 
REGRES-generated observed-minus-computed (O-C) residual (resulting in a positive correction 
to the computed observable). Hence, TSAC media calibrations will always be positive, i.e., a 
positive sign corresponds to a one-way range delay. 

2.2 Standard Mapping Tables 

Once the zenith range adjust has been accurately determined, it must be projected, or 
“mapped,” to the line of sight along the observer-to-radio source direction. In the ODP, standard 
tables of normalized scale factors are used for mapping the zenith range adjusts down to a given 
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ADJUST (DOPRNG) BY TRIG ( 31557600 . , 
0.0200, 0.0008, -0.0021, 
MODEL 

0.0870, -0.0360, 

-0.0036, -0.0002) 

-0.0336, 

0.0002, 

(WET NUPART ) FROM (72/01/01, 00 : 00) 
02:23 

DSN 

(CIO) . 

#ADJ 920121 

ADJUST (DOPRNG) BY TRIG ( 31557600 . , 
0.0033, -0.0015, 0.0005, 

MODEL 

2.0521, 0.0082, 

-0.0011, 0.0036) 

-0.0005, 

-0 . 0004, 

(DRY NUPART) FROM (72/01/01, 00: 00) 
02:23 

DSN 

(CIO) . 

#ADJ 920121 

ADJUST (DOPRNG) BY TRIG ( 31557600 . , 
0.0026, 0.0036, -0.0001, 

MODEL 

0.1149, 0.0255, 

0.0007, 0.0012) 

0.0020, 

0 . 0010, 

(WET NUPART) FROM (72/01/01, 00: 00) 
03:35 

DSN (C40 ) . 

#ADJ 920116 

ADJUST ( DOPRNG ) BY TRIG ( 31557600 . , 
0.0017, -0.0043, 0.0052, 

MODEL 

2.1579, -0.0032, 

0.0016, -0.0021) 

-0.0002, 

0.0012, 

(DRY NUPART) FROM (72/01/01, 00: 00) 
03:35 

DSN { C40 ) . 

#ADJ 920116 

ADJUST (DOPRNG) BY TRIG ( 31557600 . , 
0.0005, -0.0031, -0.0003, - 

MODEL 

0.1255, -0.0284, 

■0.0034, -0.0013) 

-0 . 0273 , 

-0 . 0094, 

(WET NUPART ) FROM (72/01/01, 00 : 00) 
03:38 

DSN ( 

C60 ) . 

#ADJ 920116 

ADJUST (DOPRNG) BY TRIG ( 31557600 . , 
0.0019, -0.0006, 0.0021, 

MODEL 

2.1094, 0.0037, 

0.0018, -0.0004) 

-0.0010, 

0 . 0036, 

(DRY NUPART ) FROM (72/01/01,00:00) 
03:38 

DSN ( C6 0 ) . 

#ADJ 920116 


Figure 3b. Format of the TSAC-supplied zenith delay tropospheric calibrations for the 
“official” seasonal model given as a 4th-degree Fourier series. 

unrefracted (straight-line) elevation angle, E, of a particular leg of the signal path between the 
observer (transmitting or receiving station) and the source (spacecraft or quasar). Separate tables 
are used for mapping the zenith dry and wet range corrections with values given at intervals of 
0.1 for elevation angles between 0 and 10 , and at intervals of 0.5° for elevation angles between 
10 and 90 . Intermediate values not given in the table are determined by the following quadratic 
interpolation scheme: 

R(E) = Rj + n (Rj+i-Rj) + 1 n (n-1) [(/? /+1 -/?,) - (R t -R, A )} (6) 

where 

R(E) = mapping factor interpolated from the standard dry or wet table 
E = Ej + A E, a given elevation angle, E , <E<E i+ 1 

n = 

Ej+\ - E, 
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The frequency of values given in the mapping tables is claimed to limit the error introduced by 
the interpolation to less than 0.1% of the total tropospheric effect [Ondrasik, 1970], 

The mapping tables currently in use were derived from ray trace results based on best-fit 
refractivity profiles from Chao’s examination of the 1967 and 1968 radiosonde balloon data 
described earlier [Chao, 1971], and normalized by the ray-traced zenith range correction. The 
ray tracing algorithm that Chao used to compute the tabular values is described in Chao [1972a], 
and in greater detail in Miller et al. [1971]. These mapping tables take into account the 
geometrical bending effect of the signal path, which becomes significant at low elevation angles 
[Chao, 1972b]. (It should be noted that these “modified” tables were first used operationally to 
support the Mariner-Venus-Mercury 1973 (MVM'73) mission and superseded the original 
mapping tables used for MM'71, which did not include the bending effect [Mottinger, 1984].) 
Although there are provisions in the ODP for overriding these tables, stored as vector arrays 
“TABDRY” and “TABWET,” the modified Chao tables have been handed down from flight 
project to flight project over the years and are still in operational use today. A standard file 
convention for the improved mapping tables was adopted in 1977 to support the Voyager 
mission (and every planetary mission since) and is provided in the Appendix. These tables are 
used by program link EDIT to compute the actual tropospheric mapping corrections. 

For use in approximate analysis, the mapping tables and quadratic interpolation algorithm 
were found to be somewhat cumbersome to implement and not readily portable to different 
computer platforms. In the development of a simple approximation model to perform the 
mapping, Chao recognized that a flat Earth model such as R(E) = 1/sin E would be inadequate, 
particularly at low elevation angles. Therefore, he extended the flat Earth approximation in a 
two-term continued fraction which could more accurately perform the dry and wet mapping: 


Ri(E) = 


sin E + 


A 

tan E + B, 


where, for i = dry: 


A dry = 0.00143 


B dry = 0.0445 


and, for i = wet: 


A wcl = 0.00035 


S we , = 0.017 


( 7 ) 


(8a) 


(8b) 


The mapping coefficients, A, and 5„ were derived from empirical fits to the original ray trace 
results. First evidence of Chao’s two-term continued fraction form of the dry mapping function 
[Eqs. (7) and (8a)] can be found in an internal document published in late 1970 [Chao, 1970]. 
The tan E term in the continued fraction was intended to be a normalization factor in order to 
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ensure that the mapping ratio would numerically approach unity at zenith, i.e., R( 90°) = 1 . Chao 
determined that the semi-empirical model was good to within 1% from the ray trace range 
correction for an elevation angle higher than 1 . It is this formula which is used to compute the 
parameter partial derivatives (partials) in the ODP [Moyer, 1984], and which is frequently 
referenced throughout the literature. Until about 10 years ago, Chao’s semi-empirical mapping 
function was arguably the simplest and most widely used closed form mapping function for any 
range of elevation angles. 

2.3 Creating the Line-of-Sight Tropospheric Correction 

A concise formula is used to determine the dry and wet tropospheric corrections (to the 
computed observable) for any radio metric data type derived from Doppler tracking, ranging, or 
VLB I techniques, assuming the zenith range adjusts and mapping factors have been adequately 
determined. This formula, described in detail by Christensen [1979a], can be expressed in 
algorithmic form as 


TROPOSPHERIC ADJUST = J TRPSGN * UC * POLY(TIME) * ELVFAC (9) 

all 

times 

(legs) 


The tropospheric sign factor TRPSGN can be ±1 depending on the direction of the leg (up or 
down) and the data type being used. The variable UC is the data type unit conversion factor 
from meters to correct units, and POLY(TIME) is any of the legitimate user-selectable 
ADJUST polynomials in meters of delay, as a function of time; these are generally the 
tropospheric zenith corrections supplied by TSAC based on the seasonal model. When utilizing 
the official seasonal model, the user should select the TRIG option (see earlier remarks). The 
variable ELVFAC is the elevation-angle-dependent mapping factor used to obtain the correction 
for the true (unrefracted) elevation angle. As described earlier, it is obtained by quadratic 
interpolation of either the TABWET or TABDRY array with the elevation angle for the 
appropriate station and time. In general, the selected model correction must be computed at 
several different times for a given data point, as the tropospheric effect can vary with time. Of 
note, the only difference between computing the model corrections for the dry versus wet 
components of the tropospheric effect are in the interpolation tables TABDRY and TABWET 
used to compute ELVFAC [ Collier , 1981], 

2.4 Artificial Deweighting of Low Elevation Data 

Although not officially part of the ODP’s tropospheric calibration model, an empirically- 
derived data editing scheme has been historically used to artificially “deweight” radio metric data 
acquired at low elevation angles [Christensen, 1979b]. The deweighting scheme was originally 
developed in an attempt to alleviate errors induced by mismodelings of the tropospheric and 
ionospheric corrections, but optimized for 2.3 GHz (S-band) data. Alternative data weighting 



algorithms for 8.4-GHz (X-band) and 32-GHz (Ka-band) Doppler data have been investigated by 
Ulvestad and Thurman [1992] and by Ulvestad [1992], Although better optimized for these 
higher link frequencies, these algorithms have yet to be incorporated into the operational ODP. 

The additional data editing process, as with the tropospheric (and ionospheric) corrections, is 
handled by the program link EDIT, and invoked with the CSP command 

WEIGHT (data-type) SIGMA (S) 

where S denotes the 1 -sigma measurement uncertainty associated with a given data point. The 
data-type argument can be any of the allowable data types supported by the ODP, e.g., radio 
metric, angle, or optical [DPTRAJ-ODP User’s Reference Manual, Vol. 1, Rev. F, 1994], Let G, 
represent the input “uncorrected” measurement uncertainty S, associated with the ith data point, 
and let C7, denote the associated “corrected” sigma. For the elevation-dependent radio metric data 
of Doppler, range, and VLBI, the weighting algorithm is given by 

WEIGHT = -V (10) 

cr, 2 

where 

Gi = Oj* E * F (11) 


Here, 


E = 1 + 05 — r 

(£»/„,+ 1 ) 

E m in, ~ minimum elevation angle (in degrees) associated with the z'th 
data point 



TCi = count time (in seconds) assigned to the z'th data point 


The count time factor, F, is to be applied to Doppler and narrowband VLBI; for range and 
wideband VLBI, F is taken to be unity. Again, variations on this weighting scheme for 
application to other elevation-dependent radio metric and non-radio metric data types are 
described in DPTRAJ-ODP User's Reference Manual, Vol. 1, Rev. F [1994]. 
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SECTION 3 


PROPOSED MODEL IMPROVEMENTS 


As seen from the previous discussion, the ODP infrastructure for computing tropospheric 
corrections is non-trivial, a result largely of early interface design requirements. Ideally, 
tropospheric and other media calibrations should be applied directly to the computed observables 
and parameter partials during the light time solution. In this way, no post-production data editing 
would be required and future modeling upgrades could easily be incorporated into the standard 
model [Collier, 1993], Unfortunately, a complete redesign of the ODP’s media calibration 
process does not appear fiscally realizable, at least not in the foreseeable future. Nevertheless, 
there are still improvements which can be incorporated into the current infrastructure, 
specifically, with respect to the fidelity of the TSAC-supplied zenith delay calibrations, and to 
the zenith delay mapping process. The proposed model improvements presented herein could be 
made with modest resource requirements and in a relatively short period of time. 

3.1 Zenith Delay Calibrations 


The TSAC software has the capability to apply real-time surface meteorological data to 
further improve the seasonal model if such a need is warranted [Runge, 1993], subject to the 
quality of the instruments used at the DSN. When surface data are used, the residuals (surface 
model minus seasonal model) can be fit with a polynomial of degree up to 9 or a Fourier series 
with degree up to 4. The real-time zenith dry delay calibration is based on a modified form of 
the dry component of the Saastamoinen [1972] model, which accounts for the site variability of 
gravity: 


Ap z (m) = 0.0022768 -ft - 


( 12 ) 


where 


f((p, H s ) = 1 - 0.00266 cos 2q> - 0.00000028 H s (13) 

Here, (p is the geodetic latitude of the DSS in degrees and H s the station mean sea level (MSL) 
height (i.e., height above the reference ellipsoid) in meters. Recall that P 0 is the surface 
barometric pressure measured in millibars. 

This is a significant improvement over the seasonal model in that it relaxes the assumption of 
constant gravitational acceleration assumed by the Berman 70 model. In addition, this model 
was derived using the refractivity constants of Thayer [1974] which supersede the obsolete 
values used by Saastamoinen and Berman. A more detailed discussion and derivation is given 
by Davis et al. [ 1 985] and Elgered [ 1 993], 
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It is important to note that Eqs. (12) and (13) model the zenith dry delay at the location of the 
barometer, therefore, an additional adjustment of the dry delay calibration must be made to the 
altitude of the local DSS. The following adjustment model performs that function and utilizes 
average surface pressure and temperature values measured at each DSCC. 

Ap Zd[y (m) = -0.0000776 I^MA//, (14) 

\ Tavg I 

This altitude adjustment model was derived from a simple application of the hydrostatic equation 
and perfect gas law [Runge, 1993]. The altitude correction, A H s , is simply the altitude offset 
between the local DSS and the reference location of the 70-m antennas at each respective DSCC. 
The major sources of zenith dry delay errors and their magnitudes are summarized below in 
Table 2 for both normal and extreme wind conditions. 


Table 2. Major sources of zenith dry delay errors under normal and extreme wind conditions. 


Error Source 

Zenith Dry Delay Errors (mm) 
Normal Extreme Wind 

Pressure sensor error 

2.3 

2.3 

Uncertainty in physical constants 

1.4 

1.4 

Use of P avg /T avg for 100 m A H s 

0.8 

0.8 

2-m error in A H s 

0.5 

0.5 

Lack of hydrostatic equilibrium 

0.2 

4.0 

Total (root-sum-square) 

2.9 

4.9 


For the real-time zenith wet delay calibration, the TSAC software utilizes the Callahan 
[1973] model, given by 


(m) = 1035 % 

io 


(15) 


where e 0 is the partial pressure of water vapor at the surface in millibars and T 0 is, as before, the 
surface temperature in kelvin. Since the TSAC software does not incorporate dew point data, it 
calculates e 0 based on a variation of the simple model of Bean and Dutton [1966]: 

eo = 6.11 (RH)o io 7 ^ 7 ®' 273 15 ^ 7 ®" 35 ' 85 ^ (16) 

where (RH) 0 is, as defined in Eq. (2), the surface relative humidity, 0 < (RH)o < 1 . The observed 
mean of e 0 for the DSN is about 8 mbar [Runge, 1993]. There is a provision to utilize the zenith 
wet delay model of Chao as an option [Chao, 1973]; however, the Callahan model is preferred 
due to its simplicity. Chao’s model for the zenith wet delay can be expressed by 
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Ap Zwc , (m) = 1 .63 + 2.05w ^ - 

To 7 ? 


( 17 ) 


where w is, as previously defined, the temperature lapse rate in K/km. The sign of w is 
determined by the linear lapse function T= To - w(h - ho), whereby the altitude, h, is constrained 
to lie between the altitude at the surface, ho, and the altitude where temperature becomes nearly 
constant as h increases (~10 km). 

It is widely known that the wet component of the troposphere is highly variable due to spatial 
and temporal inhomogeneities in the atmospheric water vapor content. Consequently, the 
Callahan and Chao models are only accurate to about 30% and correspond to an uncertainty in 
the wet delay on the order of 6 cm to 10 cm. Investigators have shown that the zenith wet delay 
model of Saastamoinen [1972] gives a modest improvement over both the Callahan and Chao 
models in certain applications (see e.g., Ifcidis [1986]). Saastamoinen’s model for the zenith wet 
delay is given by 


Ap Zwcl (m) = 0.002 2 7 7e 0 


(1255 +0.05 

To 


(18) 


A host of other zenith wet delay calibration models have been derived by a number of authors, 
including Hopfield [1971], Berman [1976] (in various forms), Ifadis [1986], and more recently, 
Baby et al. [1988], Because the variability of the atmospheric water vapor content is highly 
dependent on site and season, a universal model for the zenith wet delay has yet to be developed. 

Other more esoteric techniques include the use of Water Vapor Radiometers (WVRs), which 
infer the wet delay from sky brightness temperature measurements [Winn et al., 1976; Resch , 
1984], Theoretically, WVRs can not only accurately predict the zenith wet delay, but the entire 
wet tropospheric effect along the line of sight [Davis et al, 1985; Wilcox, 1992], Unfortunately, 
the WVR instruments have proven to be expensive and unreliable (see study by Elgered et al 
[1991]). 

The most promising approach to determining the wet component of the zenith delay appears 
to be estimating it along with the radio metric tracking data. This has been demonstrated 
successfully using space-based geodetic systems such as VLBI and the Global Positioning 
System (GPS) [Herring et al., 1990; Tralli and Lichten, 1990]. A comparison of estimates 
determined from both systems was recently published by Tralli et al. [1992], The TSAC activity 
is currently studying the feasibility of incorporating GPS-based measurements into the zenith wet 
tropospheric delay calibration process [R unge, 1994], Although not yet approved for operational 
use, the GPS technique treats the “total” zenith delay as a nuisance parameter in the data 
estimation process, and has yielded estimates of the zenith tropospheric delay to better than 1 cm 
in recent demonstrations [Lichten, 1994], The current TSAC plan is to obtain the time series of 
the GPS-derived total zenith delays over the time span of interest, and then extract the wet 
component from the total delay by simply subtracting the dry component determined from real- 
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time surface barometric pressure measurements. An actual data demonstration of this technique 
is currently being devised. 

3.2 Mapping Phase 

It is important to note that the standard mapping tables used in the ODP represent a 
“snapshot” of the refractivity profiles for a specific set of sites and seasons. Additionally, there 
is no provision in Chao’s semi-empirical mapping formula for utilizing surface weather data at a 
particular DSS site. Rigorously speaking, a new set of tables, and therefore a new “fitted” 
mapping function, should be issued for a different refractivity profile [Chao, 1972a]. 
Furthermore, Chao’s ray tracing computer programs, used to generate the ODP mapping tables, 
were re-examined by Lanyi [1983], Lanyi found that while the ray tracing algorithm was 
adequate, Chao’s assumption of the refractivity profile appeared to contradict the physical 
model, albeit at a relatively low error level, resulting in about a 2-cm discrepancy at 6° in the 
modified tables. In addition, a recent more thorough examination of the dry mapping tables 
revealed a small discrepancy of ~0.5 cm in the 1 T-elevation regime between the tabular values 
and the semi-empirical formula. It is believed that this discrepancy can be attributed to 
limitations in the original ray tracing quadrature procedure; however, no similar behavior was 
evident for the corresponding wet mapping table. 

Another important observation is that the empirical formula expressed by Eqs. (7) and (8) 
actually corresponds to the original mapping tables and not to the improved mapping tables. 
(Recall that the improved mapping tables included the geometrical bending effect.) According to 
Chao [1994], a new set of unpublished empirical constants were derived for the dry component 
of the delay around the 1977 to 1978 timeframe. The following “revised” constants enabled the 
dry approximation formula to more closely match the refined tables: 


A dry = 0.00147 


B dry = 0.0400 


(19) 


A cursory examination of the revised dry mapping function against the current dry mapping table 
does indeed suggest a marked improvement over the original formula [cf, Eqs. (7) and (8a)]. 
However, it should be noted that there are certain regions in which the revised mapping function 
exhibits a greater departure from the refined tables than the original mapping function. This 
discrepancy is most prominent at the extremely low elevation angles (< 2°), as demonstrated in 
Fig. 4, and to a lesser extent in regions between 10° and 30° (not shown in the figure). It is 
important to note that although the revised constants yield improved results over the original 
constants, the Chao approximation model is still limited by the fact that it cannot be readily 
adaptable to varying atmospheric conditions. 

As indicated in the introductory remarks, the process of mapping zenith tropospheric delay 
calibrations down to lower elevation angles has been studied extensively in recent years. Much 
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Figure 4. Differences between mapping ratios predicted by Chao’s original and revised dry 
semi-empirical approximation formulas and Chao’s improved dry mapping table 
(TABDRY) at the lower elevation angles. 

of the work in this area has come from the VLBI science community, as the refractive and delay 
effects of the neutral atmosphere effects continue to be a limiting error source in high-precision 
space-based geodetic analysis. VLBI and other interferometric techniques on long baselines 
have consistently been plagued by the fact that mutual visibility of either a spacecraft or EGRS 
between co-observing stations is limited, and therefore requires observations to be taken at 
relatively low elevation angles [MacMillan and Ma, 1994], This is certainly the case for the long 
intercontinental baselines of the DSN. The need to acquire interferometric measurements at low 
elevation angles has the advantage, however, in that such observations help reduce the 
correlation between estimates of the zenith tropospheric delay and station vertical position. This 
is arguably the principal reason that radio source observations are being acquired at ever 
decreasing elevation angles in the space-based geodetic applications. 

The remainder of this report is dedicated to assessing a set of modem and robust tropospheric 
mapping functions intended to replace the outdated Chao tables in the ODP. Although some of 
these candidate mapping functions accurately model the tropospheric mapping down to as low as 
2°, antenna pointing restrictions imposed on the DSSs due to ground obstructions, system noise, 
and signal multipath usually preclude tracking below 6°. Furthermore, spatial and temporal 
inhomogeneities in the atmosphere begin to dominate the uncertainty in the overall tropospheric 
effect in the extremely low elevation regimes (< 5°), and contribute more to the total calibration 
model uncertainty than simply the zenith delay and/or mapping errors [Lanyi, 1984], 
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SECTION 4 


MATHEMATICAL DESCRIPTION OF CANDIDATE MAPPING FUNCTIONS 


The simplest way to characterize the one-way range tropospheric correction is in terms of 
separate dry and wet components: 

AptropC^) = PZj r v^dry(£') Pz wc , ^wet(^) (20) 

Here, pz iry and pz wcl are, as before, the range corrections at zenith due to the dry and wet 
components of the troposphere, respectively. The functions /?dry(£) and ^wet(£) give the ratios 
for mapping the zenith range corrections down to the given unrefracted elevation angle, E , along 
the station-to-spacecraft (or EGRS) line of sight. Recall that in the ODP, these ratios, or 
“mapping factors,” are obtained by quadratic interpolation of Chao’s standard mapping tables, or 
by his two-term continued fraction for the parameter partials. The definition of Eq. (20) lends 
itself to easy computation of pz dry and pz wct as adjustment parameters if the user wishes to “tune” 
these parameters in the data filtering or estimation process. Hence, 

3 -£SS!^L , Rd ^E) ( 21 ) 

dpZdry 

and 

dp y (E) = «wet(f) (22) 

OPZ wcl 


Since the mapping factors, Rdry(E) and P wet (£), are functions of the incident elevation angle 
of the signal path, they can be computed by numerical ray tracing methods or can be 
approximated analytically. Ray tracing offers greater accuracy, particularly in the low elevation 
angle regimes where the tropospheric effect is difficult to model, but is often time consuming 
and impractical under varying atmospheric conditions. Often what is done in practice is to 
derive “fitted” empirically based formulas from numerical ray traces of either an idealized 
refractivity profile or from actual radiosonde profile data. Another less traditional approach is to 
derive an analytic approximation of the tropospheric mapping based on an idealized refractivity 
profile. 

Marini [1972] showed that the tropospheric mapping could be very accurately modeled by a 
continued fraction expansion of the form 


m 


sin E + ■ 


a 


sin E + 


sin E + - 


sin E + ■ ■• 


(23) 
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where the constants, a, b, c, ... , give the elevation dependence of the tropospheric delay. The 
more terms carried out in the expansion, the better the approximation of the geometrical bending 
effect. Indeed, most of the candidate mapping functions presented in this section are based on 
this formula, with the empirical constants parameterized for surface meteorological 
measurements. Consequently, these “improved” mapping functions can be adapted to a wide 
variety of applications under varying weather conditions. 


Recall that in Chao’s two-term continued fraction model [cf, Eq. (7)], the normalization 
factor, tan E, was used in order to ensure a mapping factor of unity at zenith. This is consistent 
with the simple definition given by Eq. (20). The mapping function of Eq. (23), however, does 
not ensure this property. Therefore, a ‘normalized form of the Marini mapping function is often 
used instead, which takes the form 


R(E) _ 1 + a/(l +6/(1 + c/(l + • • •))) 

sin E + ^ 

sin E + b 

sin E + £ 

sin E + ■ ■ ■ 


(24) 


This normalization process is preferred over the use of tan E since the latter can introduce a bias 
in the mapping function due its slower convergence property over a direct sin E calculation 
[Davis et al., 1985], In addition, a direct sin E calculation is not subject to an inherent analytic 
singularity at zenith (as is tan E). 


In this section, five highly accurate tropospheric mapping models are presented which are 
considered to be the state of the art. Of the five candidate models described herein, only one, 
the Lanyi mapping function [Lanyi, 1983; 1984], does not explicitly separate the dry and wet 
components, but rather combines them into a single mapping function based on an analytic 
expansion of the tropospheric delay in terms of moments of the dry and wet refractivity. The 
remaining candidate models were empirically derived and include the “CfA-2.2” of Davis et al. 
[1985], the Ifadis global solution model [Ifadis, 1986], the “MTT” functions of Herring [Herring, 
1992], and the NMF functions of Niell [Niell, 1993; 1994a; 1994b]. It is important to note that 
all of the mapping functions presented in this report (both current and proposed) were derived 
under the assumption of a spherically symmetric atmosphere, i.e., none of the mapping functions 
implicitly or explicitly account for the effects of azimuthal asymmetry. 

4.1 Lanyi Mapping Function 

In analyses of intercontinental VLBI data, it was found that the Chao semi-empirical 
mapping function was inadequate for high-precision experiments. In 1983, G. E. Lanyi 
presented an analytic tropospheric mapping function which was found to be more accurate than 
Chao s tables and analytic approximations, as well as a variety of other earlier mapping 
functions, above 4° elevation [Lanyi, 1983], These results were published later in an external 
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report [Lanyi, 1984]. Although there have been many atmospheric modeling improvements in 
recent years, Lanyi’s method arguably appears to be the most robust and complete model of its 
kind, adaptable to a wide variety of experiments involving radio interferometry and GPS-based 
analysis [McCarthy, 1992]. Its principal advantage is that the mathematical expansion used in 
the derivation is valid for any laterally homogenous atmospheric model of refractivity. Lanyi’s 
analytic approximation is estimated to be in error ( i.e ., differ from ray trace calculations) by less 
than 0.02% for elevation angles larger than 6°, less than 4 mm at 6°, and approximately 0.004% 
or 0.3 mm at 20° [Lanyi, 1984], 

In its simplest form, in which no meteorological information is available, the Lanyi mapping 
function employs standard, or “average,” values of atmospheric model parameters. Provision is 
made for specifying surface meteorological data acquired at the time of observation, that may 
override the nominal values of the model parameters. The dry and wet components of the 
tropospheric mapping are not in the Lanyi model but rather combined into a single nonlinear 
mapping function by expanding the dry, wet, and bending contributions to the atmospheric delay 
in terms of moments of refractivity. The following compact mathematical description of the 
Lanyi mapping function has been adapted from Sovers and Jacobs [1994] but is repeated and 
expanded here for completeness. 

The Lanyi one-way range tropospheric correction is defined by 


ApLanyi(£) = F(E) / sin E (25) 

where 


F(E) = Pz^dryiE) + Pz^ME) 

+ [pZ«yF bendlCE) + 2pZ d , y Pz wc / 7 bend2(£) + PzL Ebcndl(E)} I A (26) 

+ Pi^bend4(£) I A 

The quantities p Zdry and pz wct are the usual zenith dry and wet tropospheric delays; A is the 
atmospheric scale height given by A = kTJmg c , with k the Boltzmann’s constant in erg/K, To the 
average surface temperature in kelvin, m the mean molecular mass of dry air in grams, and g c the 
gravitational acceleration at the center of gravity of the air column in erg/g’cm. Using the 
standard values k = 1.38066 x 10' 16 erg/K, m = 4.8097 x 10' 23 g, g c = 978.37 erg/g*cm, and the 
average temperature for DSN stations To = 292 K, the atmospheric scale height is computed to be 
A = 8567 m. 

The dry, wet, and bending terms are expressed as of moments of the refractivity: 

F dry (E) = A io(£) GiXM, io, u ) + 3 guM 2W G i 0 , u ) / 2 (27a) 

Fwei{F) = Ao\(E)G\XMm! Mw\, u) / A/ooi (27b) 
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^bendl(£) = \oG\Mno, w)/sin 2 E - M 02 oG 3 (A/i 2 o/Mo 2 o, w )] / 2 tan 2 E (27c) 

^bend2(£) = -MonG \M\uIMqu, u ) / 2 Moot tan 2 E (27d) 

T'bendsC^) = -A/ 002 G (M 102 /A/ 002 , w ) / 2 Mooi tan 2 £ (27e) 

^beiKMCE') = M 030 G 3 (Mi 3 o/M) 3 o, U )/ tan 4 £ (27f) 


Note that in Za/iyi [1984], there are two misprints in the corresponding equations. In Eq. (3) of 
the referenced document, the misprinted factor of 3/4 should be replaced by 3/2, while in Eq. (5) 
of the referenced document, the erroneous multiplier of -1/2 in the Fb en d4(£) term should be 
removed [Lanyi, 1993], These changes are reflected in the above equations. G(q,u ) is a 
geometric factor given by 


with 


G(q,u) = (1 + qu) m 
u = 2<t/ tan 2 E 


(28) 

(29) 


where cr= A/R is a measure of the curvature of the Earth’s surface with standard value 0.001345. 
An approximation for the local radius of curvature, R, is given in Eq. (37) below. 


The quantities A/ m (E) and Mq m are related to moments of the atmospheric refractivity. A I0 (E) 
involves the dry refractivity, while A 0 i(E) is the corresponding wet quantity. The A/ m (E) are 
given by 


10 


n=l k=0 


(- 1 ) n+k (2n- 1 ) ! \M n .kj<m 

u 

\XM\im 

2 n k\(n-k)\ 

_ 1 + XuM U m /MQim _ 

1 

s 

I 


(30) 


with the scale factor A = 3 for E < 10° and X- 1 for E> 10°. Only the two combinations 
(l,m) = (0,1) and (1,0) are needed for the Ai m (E). The moments of the dry and wet refractivities 
are defined as 


M nij = J dq <7"/dry(?)/wet(?) (31) 

Jo 

where _/dry >wet (gO are the surface-normalized refractivities. Here, n ranges from 0 to 10, i from 0 
to 3, and j from 0 to 2; not all combinations are needed. Carrying out the integration in Eq. (31) 
for a three-section temperature profile gives an expression for the general moment, M ni j : 
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M, 


ti 

r m y/«! = (1 - e- fl 9')/o" +1 + e- a *' [l - Tt n+ \quqi)] ]1 7"? 


7= o b+i + 1 


+ e -^,^ + « + i( 9l)92 )/ a «+i 


( 32 ) 


Here, 

T 2 {q\,qi) = 1 - (^2 - <?i) / a (33) 

The quantities gi and g 2 are the scale-height normalized inversion and tropopause altitudes, 
respectively. For the standard atmospheric model, q\ = 0.1459 and q 2 = 1.424. The constants a 
and b are functions of the dry ( a = 5.0) and wet (/J = 3.5) model parameters, as well as of the 
powers of the refractivities (/ and j) in the moment definitions. Table 3 gives the necessary a's 
and b ' s for the current mechanization (version “42map” 30-sep-85). Note that the normalization 
is such that Af 0 io = 1; this moment has therefore been excluded from Eqs. (27a) through (27f). 


Table 3. Dependence of constants a and b on the Lanyi tropospheric model parameters a and 0 

(version “42map” 30-Sep-85). 


i 

j 

Moment 

Type 

a 

b 

i 

0 

dry 

i 

a - 1 

0 

1 

wet 

p 

a(5 - 2 

2 

0 

dry squared 

2 

2(a- 1) 

1 

1 

dry times wet 

0+1 

(5(a+ 1)-3 

0 

2 

wet squared 

2(5 

2(a(5 - 2) 

3 

0 

dry cubed 

3 

3(«-l) 


Provision is made for input of four meteorological parameters to override the nominal 
(average) values of the Lanyi model. These are: 1) the surface temperature, T 0 , which 
determines the atmospheric scale height (nominal value 292 K); 2) the temperature lapse rate, w, 
which determines the dry model parameter a (nominal values w = -6.8165 K/km, a = 5.0); 3) the 
inversion altitude, hi, which determines q\ = h\/A (nominal value h , = 1.25 km); and 4) the 
tropopause altitude, h 2 , which determines q 2 = h 2 / A (nominal value h 2 = 12.2 km). A fifth 
parameter, the surface pressure, P 0 , is not currently used. Approximate sensitivity of the 
tropospheric delay (at 6° elevation) to the meteorological parameters is -0.7 cm/K for surface 
temperature, -2 cm/(K/km) for lapse rate, and -2 cm/km for inversion and 0.5 cm/km for 
tropopause altitudes, respectively. 

In addition, there are two site-dependent parameters which may be input to override the 
standard values. The first is the gravitational acceleration at the center of gravity of the air 
column, g c , defined earlier with a nominal value of 978.37 erg/g*cm. Recall that it is used, along 
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with surface temperature, To, to determine the atmospheric scale height A. The gravity term, g c , 
can be approximated as [ Lanyi , 1983]: 

gc ~ go - 2(kTo/mR)Mo\o (34) 

The second temi in the difference represents the “shift” in gravity (-2.25 erg/g*cm), which is a 
function of surface temperature and temperature profile parameters. The approximate surface 
value ofg (erg/g*cm) comes from Saastamoinen [1972]: 

go = 980.62 ( 1 - 0.00265 cos 2 (p - 0.0003 1 H c ) (35) 

where (p is the geodetic latitude of the station in degrees, and H c is the height in kilometers of 
the center of mass of the vertical column of air. The default value forgo is taken to be 980.62 
erg/g*cm, corresponding to gravitational acceleration at sea level ( H c = 0) and midlatitude 
(<jP=45°). Because of the difficulty in implementing Eq. (34), Saastamoinen’s “average” 
approximation formula is often used for computing the site-dependent variability of gravity [cf, 
Eq- (13)]: " 


g c = 978.4 (1 - 0.00266 cos 2 (p - 0.00028 H s ) (36) 

where H s is the MSL height of the station in kilometers. 

The other site-dependent quantity, which is used to determine the variation of go as well as 
the Earth curvature parameter, a, is the local radius of curvature, R, which can be approximated 
as 

R « R e (1 - 0.003353 cos 2<p+ 0.000157 H s ) (37) 

Here, R e is the mean equatorial radius of Earth (default value 6378.16 km). 

In analysis of data for which meteorological parameters are not available, there is an 
approximate form of the mapping function (again, version “42map”) which involves a one- 
parameter estimate. This parameter, simply denoted here as p, accounts for deviations from 
standard meteorological conditions. It is a modified form of the ~l/tan 2 E term but includes wet 
a priori shift values and is more realistic at low elevation angles than 1/tan 2 E. The bending 
effect is not included in this approximation, as it is assumed to be modeled sufficiently by 
solving for the zenith dry plus wet delay, hence, the dry and wet atmospheric refractivity factors 
Aio(E) and Aoi(E) are taken to be unity. For this special case, the one-parameter correction is a 
fraction of the approximate tropospheric correction to (pz iry + Pz wct ) / sin E: 


A/>L,n,i(0 = <p 4 „ + PZJ / ™ E + p - APL " yl<£) 

dp 


(38) 
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The partial derivative term is given by 


^ApLanyiC^O 

dp 


(PZdrv + PZwJ^'IO 

G(Mno, «)[1 + G(Mno, «)]sin E 

Pz wcl w(Afuo - Mm/Mood 

G(Muo, u)G(Mm/M m , w)[G(A/iio, u) + G(Mm/M m , w)]sin E 


(39) 


Because the definition of the one-way range tropospheric correction for the Lanyi model does 
not explicitly separate the dry and wet components, as in Eq. (20), the parameter partials for the 
conventional tuning parameters, p Zin and pz wel , are expressed as 


— p Lanyi (£) = [ Fdry(£) + 2 p Zir f hmdl (E) / A ] / sin E 

^PZdry 


+ 


[2p Zw Aend2(£) / A + 3 p z ; y F ben d 4 (E) / A J / sin E 


(40) 


y A EL = [F wet (E) + 2p Zi!y F bend2 (E) / A + 2p z ^F hcnd3 (E) / A ] / sin E 
dpz~„a 


(41) 


4.2 Davis et al. (“CfA-2.2”) Mapping Function 

Around the same time that Lanyi published his new mapping function, scientists from the 
U.S. East Coast VLBI community were developing a new empirically based formula for the dry 
component of the tropospheric mapping. The resulting mapping function, designated “CfA-2.2,” 
was introduced by J. L. Davis and his colleagues in 1985 and resulted from extending Chao’s 
empirical formula in the form of a three-term Marini continued fraction, fit to ray traces of an 
idealized atmospheric model of refractivity [ Davis et al., 1985]. Additionally, the mapping 
function was parameterized for surface meteorological measurements; namely, surface pressure 
and temperature, relative humidity, lapse rate, and tropopause altitude. A companion mapping 
function for the wet component of the delay was not derived. The CfA-2.2 mapping function is 
estimated to deviate from ray trace calculations by less than 5 mm for elevation angles down to 
5° [Davis et al., 1985]. 

The mapping function, denoted here as Rq fA-2.2(£)> is given by 

*CfA-2.2(£) = L — (42) 

sin E + 

tan E + 

sin E + c 
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where the three parameters a , b, and c are expressed in terms of meteorological data as: 


a = 0.001 185 [l + 0.6071 x 10' 4 (ft, - 1000) 

-0.1471 x 10' 3 e„ + 0.3072 x 10' 2 (ft -293.15) (43a) 

+ 0.1965 x 10' 1 (w + 6.5)- 0.5645 x 10' 2 [hi- 11.231)] 

b = 0.001 144 [l +0.1 164 x 10' 4 (ft, - 1000) 

+ 0.2795 x 10' 3 eo + 0.3 109 x lO' 2 (ft -293.15) (43b) 

+ 0.3038 x 10' 1 (w + 6.5)- 0.1217 x 10 1 (h 2 - 11.231)] 

c = -0.0090 (43c) 

Here, ft, is the total surface pressure measured in millibars, and e<, is the partial pressure of water 
vapor at the surface, also measured in millibars. Nominal values are taken to be Po = 1000 mbar 
and eo = 0. ft is the surface temperature, originally given in degrees Celsius but shown here in 
kelvin for consistency with the Lanyi mapping function, w is the temperature lapse rate (by 
definition, vv < 0) in K/km, and h 2 is the altitude of the tropopause in kilometers. The nominal 
values for these parameters are ft = 293.15 K (- 20 °C), w = -6.5 K/km, and h 2 = 1 1.231 km. 
Approximate sensitivity of the tropospheric delay (at 5° elevation) to the meteorological 
parameters is [Davis et al., 1985]: -0.17 cm/mbar for surface pressure, -0.75 cm/K for surface 
temperature, -4.4 cm/(K/km) for lapse rate, and 1.1 cm/km for tropopause altitude. 

When evaluating the performance of this new mapping function with VLBI data for a 
multistation intercontinental geodetic experiment, Davis et al. used the CfA-2.2 to represent both 
the dry and wet mapping components of the delay. By not using a mapping function derived 
explicitly for the wet component, a small error was introduced in the VLBI results, as the wet 
zenith correction was multiplied by the difference between the CfA-2.2 and the “true’’ wet 
mapping function [Davis et al., 1985]. An alternative approach sometimes used in practice is to 
utilize the CfA-2.2 mapping function for the dry component of the delay and the Chao wet 
mapping function for the wet component (see e.g., MacMillan and Ma [1994]). 

4.3 Ifadis Global Mapping Functions 

An extensive study of newly proposed global-, site-, and climate-dependent tropospheric 
calibration models was conducted by I. 1. Ifadis and published in 1986 [Ifadis, 1986]. Ifadis’ 
new models for the zenith delay calibration and tropospheric mapping were based on ray trace 
analyses of radiosonde profiles collected from various sites distributed over large areas of the 
world and under different climactic conditions. More specifically, the data base used in his study 
consisted of height profiles of pressure, temperature, and relative humidity, as well as wind speed 
and direction, which were collected over a 3 year period from 47 globally distributed sites. 
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Ifadis derived dry and wet mapping functions based on the “unnormalized” continued 
fraction form of Marini [Eq. (23)], truncated at three terms and at four terms in the expansion. 
Different mapping coefficients were derived depending on whether the global-, site-, or climate- 
dependent model was being used, but all are linear functions of surface temperature, pressure, 
and relative humidity. Because a great deal of care has to be exercised when using the climate- 
dependent model or in “optimizing” the mapping coefficients for the site-dependent model, only 
the global-solution model is referred to in this report. The Ifadis three-term global mapping 
function is estimated to be in error by less than 4 mm at 2° elevation [Ifadis, 1986], 

In order to be consistent with the definition of Eq. (20), a normalized form of the Ifadis 
global mapping function is assumed herein, which takes the form 


W£> - — - a,/(1+v ‘ 1+c ')) 

sin E + l - 

sin E + — — ^ 

sin E + Cj 


(44) 


where, for / = dry: 


a dry = 0.001237 + 0.1316 x 10' 6 (P 0 - 1000) 

+ 0.1378 x lO’^ro- 288.15) (45a) 

+ 0.8057 x 1 O' 5 VftT 

b dry = 0.003333 + 0.1946 x lO'Vo - 1000) 

+ 0.1040 x 10' 5 (r o - 288.15) (45b) 

+ 0.1747 x 1 O' 4 


c dry = 0.078 


(45c) 


and, for i — wet: 


u we , = 0.0005236 + 0.2471 x 10' 6 (P 0 - 1000) 

- 0.1724 x 10' 6 (7o - 288.15) (46a) 

+ 0.1328 x 10‘ 4 V^ 

6 we t = 0.001705 + 0.7384 x 10‘ 6 (P 0 - 1000) 

+ 0.3767 x lO^To - 288.15) (46b) 

+ 0.2147 x 10‘ 4 Veo 

£\vet ~ 0.05917 (46c) 
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A misprint in Eq. (6.10b) of Ifadis’ original report for the b diy coefficient has been corrected here 
[ MacMillan , 1994]. The meteorological parameters used in the mapping coefficients are defined 
as Po for surface pressure in millibars, To for surface temperature (shown here in kelvin) and eo 
for the partial pressure of water vapor in millibars. Nominal values for these parameters are 
Po = 1000 mbar, To = 288.15 K (= 15 °C),ande o = 0. 

4.4 Herring (“MTT”) Mapping Functions 

The “MTT” dry and wet mapping functions were introduced by T. A. Herring (co-author of 
the CfA-2.2 dry mapping function) in 1992 [Herring, 1992], The MTT mapping functions 
provide empirically based formulas for both the dry and wet components of the tropospheric 
effect and take the form of three-tenn nonnalized Marini continued fractions. The coefficients of 
the mapping function are linear functions of surface temperature (which is the only 
parameterized meteorological input) and station geodetic latitude and height. The dry and wet 
mapping functions were fit to ray traces of 2 years of atmospheric profiles from 1 1 North 
American sites. The MTT mapping function is estimated to be in error by less than 1 mm at 5° 
elevation, if the coefficients in the mapping function are known, and ~30 mm if the coefficients 
are determined from their correlation with location and surface temperature [Herring, 1992], 

The function corresponding to the dry and wet mapping is given by 




1 + dj / ( 1 + bj/( 1 + Cj)) 

sin E + — 

sin E + — 

sin E + c, 


(47) 


where, for / = dry: 

a dry = [ 1.2320 + 0.0139 cos (p - 0.0209 H s + 0.00215 (7o - 283.15)] • 10' 3 (48a) 

bdry = [3.1612 -0.1600 cos <p- 0.0331 //, + 0.00206 (r 0 - 283.15)] • 10' 3 (48b) 

c dr y = [71.244-4.293 cos <p- 0.149//, - 0.0021 (T 0 - 283.15)] • 10‘ 3 (48c) 

and, for i = wet: 

a wet - [0.583 -0.011 cos (p- 0.052//, + 0.0014 (ft -283. 15)]- 10' 3 (49a) 

fcwet = [1.402 - 0.102 cos (p- 0.101 H s + 0.0020 (To -283.15)] • 10" 3 (49b) 

c we , = [45.85 - 1.91 cos <p- 1.29//, + 0.015 (To -283.15)] • 10' 3 (49c) 
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Here, <p and H s denote the geodetic latitude and MSL height of the station, respectively, with (p 
specified in degrees and H s in kilometers; To is the surface temperature (again, shown in kelvin) 
with nominal value To - 283.15 K (= 10 °C) . 

4.5 Niell (“NMF”) Mapping Functions 

All of the candidate mapping functions discussed thus far have either been based on 
parameterizations of surface meteorology or have some direct dependence on surface 
meteorological data. A somewhat different approach was suggested recently by A. E. Niell 
[Niell, 1993; 1994a; 1994b], Niell claims that mapping functions which depend on surface 
temperature are limited in their accuracy since the variability in temperature near the Earth’s 
surface is much greater than at higher (> 2000 m) altitudes. Therefore, he based his mapping 
functions (dry and wet) on temporal fluctuations in the bulk of the atmosphere rather than on 
meteorological conditions near the surface. The Niell global mapping functions, designated 
“NMF,” were derived from temperature and relative humidity profiles of the U.S. Standard 
Atmosphere at North latitude regions 15° (tropical), 30° (subtropical), 45° (midlatitude), 60° 
(subarctic), and 75° (subarctic) for the months of January (Winter) and July (Summer). Niell 
assumes that the Southern and Northern hemispheres are anti-symmetric in time, i.e., that the 
seasonal behavior is the same. In addition, he assumes that the equatorial region is described by 
the 15° N latitude profile while the polar regions are described by the 75° N latitude profile. The 
NMF mapping functions are estimated to be in error by less than 4 mm from 12° down to 3°, 
comparable to the MTT mapping functions of Herring, but with smaller biases relative to ray 
traces than the MTT functions [Niell, 1994a]. 

The NMF mapping functions utilize the same normalized continued fraction expansion of 
Marini that was used in the Ifadis and MTT dry and wet mapping functions, truncated at three 
terms, but with considerably different parameterizations of the mapping coefficients. In addition, 
the NMF dry mapping function contains a correction term to account for a geoidal station height 
dependence, reflecting the fact that the ratio of the atmosphere “thickness” to the radius of 
curvature decreases with height [Niell, 1994b], 

The NMF dry mapping function can be expressed as 

^NMF dry (£’) = nidryiE) + AmiryiE) (50) 


where 


m dry(E) 


1 "E Odry/(l "E bfc y/(l E Cdry)) 


sin E + 


a dry 


sin E + 


bd, 


ry 


sin E + Cdry 


(51) 
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The mapping coefficients, a dxy , b dry , and c dry , are linear fixnctions of tabular latitude ty, at time t 
from January 0.0 (in UT days): 


fl dry (<Pi,t) = a dry ayg ((pi) + a dT y amp (<p,)C0S 


r-To 

.365.25 


(52) 


Similar equations exist for the b dry and c dry mapping coefficients. Here, T 0 denotes the “phase” 
in UT days past January 0.0 and not the surface temperature as in earlier mapping function 
descriptions. A nominal value of T 0 = 28 days is used in the current version of the dry 
(hydrostatic) mapping function (designated “nmfh2” 26-Jan-94), corresponding to the day that 
exhibited the maximum value of the mapping function based on data analysis. The latitude 
dependence reflects the five different latitude regions discussed earlier (at 15°, 30°, 45°, 60°, and 
75° N), thus, i ranges from 1 to 5 in Eq. (52). The average ( avg ) and amplitude {amp) 
coefficients are given in Table 4a. Intermediate values for the mapping coefficients are obtained 
by linear interpolation. 


Table 4a. Mapping function coefficients for the dry (hydrostatic) NMF mapping function 

(version “nmfh2” 26-Jan-94). 


Coefficient 



Latitude 




15” 

30” 

45” 

60° 

75” 




average (x 10 -3 ) 



a dry uvg 

1.2769934 

1.2683230 

1.2465397 

1.2196049 

1.2045996 


2.9153695 

2.9152299 

2.9288445 

2.9022565 

2.9024912 

c dry aV f> 

62.610505 

62.837393 

63.721774 

63.824265 

64.258455 




amplitude (x 10' 5 ) 


a dry amp 

0.0 

1.2709626 

2.6523662 

3.4000452 

4.1202191 

^dry amp 

0.0 

2.1414979 

3.0160779 

7.2562722 

11.723375 

^*dry am p 

0.0 

9.0128400 

4.3497037 

84.795348 

170.37206 


The correction term, A m dry (E ) , is given by 

A m^E) = (53) 

ah 

where H s is the MSL height of the station. The analytic height correction coefficient, 
dm diy (E) / dh, is taken to be 


dm dry {E?) j 

dh sin E 


f {E, a h(,b ht,c hi) 


(54) 
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Here, f(E;ah t ,bht,Cht) represents the three-term continued fraction expressed by Eq. (51) using 
the following values for the “fitted” height parameters: 


a ht = 2.53 x 10' 5 

(55a) 

b ht = 5.49 x 10* 3 

(55b) 

Cht = 1.14 x 10' 3 

(55c) 


Because only the latitude dependence is predictable for the corresponding wet mapping 
function, it takes a simpler form: 


^NMF wcl (E) = 


sin E + - 


sin E + bwet 


sin E + c, 


(56) 


The values for the a wet , b wet , and c wet coefficients are provided in Table 4b, corresponding to the 
five different latitude regions. As with the dry mapping coefficients, intermediate values not 
shown in the table are determined by linear interpolation. The current version of the wet 
mapping function is designated “nmfw2.” 


Table 4b. Mapping function coefficients for the wet NMF mapping function 

(version “nmfw2” 26-Jan-94). 


Coefficient 



Latitude 




15° 

30° 

45° 

60° 

75° 




(x 10- 4 ) 



^wet 

5.8021897 

5.6794847 

5.8118019 

5.9727542 

6.1641693 




(x lO- 3 ) 



^wet 

1.4275268 

1.5138625 

1.4572752 

1.5007428 

1.7599082 




(x 10- 2 ) 



C\vet 

4.3472961 

4.6729510 

4.3908931 

4.4626982 

5.4736038 
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SECTION 5 


EXPERIMENTAL EVALUATION 


Because ray tracing through an extensive radiosonde data set was recently conducted [see 
Mendes and Langley, 1994], no attempt was made here to perform additional ray tracing 
computations to evaluate the performance of the candidate tropospheric mapping functions 
presented thus far. The Mendes-Langley study culminated in a comprehensive analysis covering 
different climatic regions of a number of tropospheric mapping functions currently used in 
academia and industry, including all of the mapping functions described in this report. In an 
effort to be complete, this section describes a recent statistical evaluation of the candidate 
tropospheric mapping functions based on actual VLBI data acquired by the DSN. 

5.1 Overview 

Tests were devised to evaluate the performance of various tropospheric mapping functions 
using VLBI measurements with DSN antennas carried out during the past 2 decades [ Sovers and 
Lanyi, 1994], These are typically 24-hour observing sessions using two antennas separated by 
either of the two intercontinental baselines (-8,400 or 10,600 km), looking over as large a region 
of the sky as permitted by mutual source visibility. Because the extreme baseline length is a 
large fraction of an Earth diameter, the region of the sky accessible at either site is severely 
limited, and numerous observations are made at the lower elevation limits of the antennas (~6°). 
These low-elevation observations help to decorrelate zenith delays at the two stations in data 
analyses, but require high accuracy in the tropospheric mapping at low elevation angles. While 
tropospheric delay is also an important part of the theoretical model in orbit determination with 
the ODP software set [Mottinger, 1982], similar tests for the Doppler and range observables in 
spacecraft tracking are precluded by the lack of accepted observations below 15° elevation, 
which is the region most sensitive to details of the mapping algorithms. 

In order to avoid the possibility of drawing false conclusions induced by the much higher 
system noise level of the early observations (Mark II —10 cm, Mark III -2 cm), only a subset of 
the VLBI data was selected for detailed analyses. It included all the newer data, and was 
recorded with Mark III data acquisition systems. The data spanned the time period 1988 to 1993; 
a total of 1 1 ,897 delay and delay rate pairs were used. A standard VLBI parameter estimation fit 
was performed using JPL’s MODel ESTimation (MODEST) software set [Sovers and Jacobs, 
1994]. The estimated parameters included positional coordinates of 283 radio sources, a pair of 
nutation angle offsets for each session (longitude and obliquity), and station coordinates for each 
overseas station (Australia and Spain) for each of the 61 sessions. The right ascension of one 
source, the nutation model on one day, and the Goldstone station coordinates were kept fixed. A 
new value of the zenith tropospheric value was estimated every 2 hours at each station. 
Correlations among the delay and delay rate observables due to tropospheric fluctuations were 
ignored here, but will be considered in a future study. More details of both the DSN data 
acquisition and parameter estimation procedures have been previously published [Sovers et at. 
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1988; Sovers, 1991]. MacMillan and Ma [1994] recently performed a similar evaluation for the 
Chao, CfA-2.2, Ifadis, and MTT tropospheric mapping functions, using the NASA Crustal 
Dynamics and IRIS databases. 

A post-fit analysis estimates linear time variation of the station coordinates, as well as the 
Goldstone-overseas baseline lengths. Such fits and post-fit baseline analyses were repeated for 
each of ten tropospheric mapping functions, with all other aspects of modeling and parameter 
estimation (the latter with one exception) being identical. The mapping functions are listed in 
Table 5. For those functions requiring surface meteorological measurements (No.s 4 through 9) 
data were taken from the DSCC Media Calibration Subsystem database [Runge, 1993]. 
Temperature lapse rates had seasonal variations at the Australia and Spain stations [Smith, 1992], 
while the default w = -6.8165 K/km was used for Goldstone. The two remaining atmospheric 
parameters (inversion and tropopause heights) were assigned default values of h\ = 1.25 and h 2 = 
12.2 km, respectively. The “Lanyi estimated” function requires some additional explanation. It 
uses the Lanyi [1984] function to accurately map the zenith delays, based on the available 
tabulated surface and balloon meteorological data. A crude single-parameter approximation is 
then used to account for the effect of deviations between the real and tabulated meteorological 
data on the mapping function [cf. Section 4.1], One such parameter is estimated at each station 
for every observing session. 


Table 5. Tropospheric mapping functions used in fits to DSN VLB! data. 


Mapping 

Function 

JPL VLBI 
Archive ID 

Reference 

Comments 

1 = Chao original 

152 

[Chao, 1974] 


2 = Chao revised 

153 

— 

Revised constants 

3 = Chao tables 

154 

[Chao, 1977] 

Tables used in ODP 

4 = Lanyi standard 

151 

[Lanyi, 1984] 


5 = Lanyi updated 

155 

[Lanyi, 1983] 

Geometry and gravity 
curvature corrections 

6 = Lanyi estimated 

156 

[Sovers and Jacobs, 1994] 

One parameter/station 
estim. per 24-h expt. 

7 = CfA-2.2 

158 

[Davis et al., 1985] 


8 = Ifadis 

157 

[Ifadis, 1986] 


9 = MTT 

159 

[Herring, 1992] 


10 = NMF 

160 

[Niell, 1993; 1994a; 1994b] 



Examination of tropospheric mapping functions includes: 1) comparison of the functional 
forms of the various mapping functions for the DSN sites for the particular subset of 
meteorological conditions prevailing at the times of the VLBI measurements, and 2) comparison 
of VLB I delay and delay rate residuals and baseline length scatter resulting from multi-parameter 
estimation in model fits to the VLBI observables. The first category thus reflects the effect on 
the mapping functions of the elevation distributions and weather conditions at the times of VLBI 
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experiments, while the second category is a quantitative assessment of how well each mapping 
function represented the data. These two categories of tests are described in greater detail below. 

5.2 Statistics of Mapping Function Delay Values 

Direct comparisons yield voluminous data, most of which will not be presented here. They 
are in essential agreement with the results presented by Mendes and Langley [1994], 
Discrepancies between mapping functions generally increase rapidly with decreasing elevation 
angle, in most cases reaching ~10 cm of tropospheric delay at the DSN lower observational limit 
of 6 . A histogram of the elevation angles of the DSN VLBI observations is shown in Figure 5. 
The distribution is seen to peak between 10° and 15°, with more than half of the observations 
being below 30° elevation at one or both ends of the baseline. This distribution of observations 
with elevation angle, along with the magnitude of the differences between the mapping 
functions, is generally consistent with the post-fit analyses. 
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Figure 5. Histogram of elevation angle distributions of DSN Mark III VLBI observations. 


The NMF function appears to show a performance nearly equivalent to that of the Lanyi 
mapping function, therefore, greater comparative detail is given here. Figure 6 shows the NMF 
minus Lanyi standard differences as a function of elevation angle. The two functions were 
evaluated at each station for all 11,897 observations. The results were placed in 0.5° bins for 
elevations lower than 20°, and 1° bins above 20°; the error bars in Fig. 6 correspond to the scatter 
from the average in each bin. The differences are as large as 3 cm just above 6° elevation, and 
decrease toward zero when approaching zenith. 
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Figure 6. Tropospheric delay difference vs. elevation angle for DSN Mark III VLBI 
observations, NMF minus Lanyi mapping functions. (Each error bar is the 
standard deviation from the mean of all values in each elevation bin.) 


5.3 Statistics of Post-fit Residuals 

The post-fit quantities that were examined included the delay and delay rate residuals, the 
“baseline scatter” and chi square for a time-linear fit to baselines (estimated independently for 
each observing session), and the baseline length bias (change in baseline length at a given epoch 
induced by a change in the mapping function). Criteria for a good mapping function are small 
values for the residuals and the baseline scatter. Differences in the average angular dependence 
of the mapping functions, as well as in the resulting baseline scales, can point to particular 
elevation ranges or baselines which are not properly modeled. 

Table 6 shows the RMS delay and delay rate residuals and total chi squares for the 10 VLBI 
fits. Normalized £ 2 values are not shown because of the uncertainty in partitioning parameters 
between the two observable types. The “improvement” columns in this table (as well as in 
Tables 7b and 8) contain the signed RMS difference between the residuals given by each 
mapping function and Lanyi standard. Focusing first on the residual results in Table 6, the 
mapping functions are seen to fall into three broad groups. Both Chao functions and the Chao 
tables produce fits to VLBI data that are definitely inferior, especially for delay rates. All other 
functions do not differ by more than a few picoseconds in delay residuals from the Lanyi 
function and its variants. Only the NMF and Ifadis functions, and Lanyi with parameters 
estimated, improve both the delay and delay rate residuals. As may be seen additionally from the 
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Table 6. Mark III VLBI residuals. 


Mapping 

Function 

Delay (r) 
RMS (ps) X 

Delay rate (t) 
RMS (fs/s) X 2 

Improvement 
Ar(ps) Aiffs/s) 

1 

71.6655 

9665 

131.512 

13140 

-19.4 

-38.2 

2 

69.2402 

8932 

126.494 

11883 

-5.9 

-12.7 

3 

70.4522 

9176 

154.737 

21456 

-14.3 

-90.0 

4 

68.9867 

8835 

125.856 

11701 



5 

69.0001 

8838 

125.869 

11702 

-1.4 

-1.8 

6 

68.8222 

8784 

125.314 

11599 

+4.8 

+ 12.7 

7 

69.0492 

8844 

126.315 

11788 

-2.9 

-10.8 

8 

68.9374 

8837 

125.762 

11717 

+2.6 

+4.9 

9 

69.0331 

8838 

125.767 

11695 

-2.5 

+4.7 

10 

68.9273 

8813 

125.608 

11682 

+2.9 

+9.4 


subsequent results, the NMF and Ifadis functions show very similar behavior in all categories, 
and may be considered to be basically identical. 

In order to assess the behavior of each tropospheric mapping function in various elevation 
angle ranges, the delay residuals were divided into six elevation bins in Tables 7a and 7b: below 
10 , 5 bins up to 30°, and observations from 30° to zenith. Table 7a gives the raw results, while 
Table 7b relates them to the Lanyi standard function. It is seen that the NMF and Ifadis (and 
Lanyi with estimated parameters) improvement in residuals is not uniform across the range of 
elevations; they do well at very low angles, but not in the 10° to 15° and 25° to 30° ranges (this 
holds true for most of the newer functions). No extraordinary elevation-partitioned results are 
seen, with one exception: it can be noted that the Lanyi parameters-estimated function apparently 
achieves reductions in residuals by improving the mapping function “shape” at very low 
elevations. 


Table 7a. Delay residuals (ps) by elevation angle. 


Mapping 

Function 

0-10 

Range of Elevation Angles, deg 
10-15 15-20 20-25 25-30 

30-90 

1 

98.993 

76.114 

66.665 

64.213 

65.090 

66.386 

2 

94.562 

73.217 

65.754 

62.232 

62.562 

64.346 

3 

97.358 

74.320 

66.923 

62.694 

62.537 

65.477 

4 

94.117 

72.692 

65.883 

62.025 

61.914 

64.221 

5 

94.134 

72.708 

65.904 

62.036 

61.911 

64.235 

6 

93.679 

72.807 

65.770 

61.903 

61.938 

64.000 

7 

94.419 

72.716 

65.873 

61.990 

61.900 

64.267 

8 

94.049 

72.782 

65.632 

61.965 

62.171 

64.121 

9 

94.174 

72.918 

65.777 

62.027 

62.173 

64.210 

10 

93.994 

72.772 

65.741 

61.938 

62.151 

64.104 
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Table 7b. Improvement of delay residuals by elevation range. 


Mapping 

Function 

0-10 

Improvement (ps) over Lanyi standard 
10-15 15-20 20-25 

in quadrature 
25-30 30-90 

1 

-30.7 

-22.6 

-10.2 

-16.6 

-20.1 

-16.8 

2 

-9.2 

-8.8 

+4.1 

-5.1 

-9.0 

-4.0 

3 

-24.9 

-15.5 

-11.8 

-9.1 

-8.8 

-12.8 

4 

5 

-1.8 

-1.5 

-1.7 

-1.2 

+0.6 

-1.3 

6 

+9.1 

-4.1 

+3.9 

+3.9 

-1.7 

+5.3 

7 

- 7.5 

-1.9 

+1.1 

+2.1 

+1.3 

-2.4 

8 

+3.6 

-3.6 

+5.7 

+2.7 

-5.6 

+3.6 

9 

-3.3 

- 5.7 

+3.7 

-0.5 

-5.7 

+ 1.2 

10 

+4.8 

-3.4 

+4.3 

+3.3 

-5.4 

+3.9 


Table 8 shows baseline length results from the VLBI data analyses. It is seen that, without 
exception, any departure from the standard Lanyi mapping function increases the baseline scatter 
and x 1 P er degree of freedom (xl) (with the exception of CfA-2.2 on the 15-65 baseline); 
sometimes by substantial amounts. The statistical significance of these differences can be 
inferred from the formal uncertainties of xl ■ The number of degrees of freedom N in the 
Califomia-Australia (15-45) and Califomia-Spain (15-65) fits are 30 and 27, respectively. 
Thus c x i » (2!N) m = 0.26 and 0.27. Even for the best fits, however, xl shows substantial 
departures from unity. This originates both from model inadequacies and from underestimated 
observable errors. It is assumed that the best fits correspond to the best modeling, and that in this 
case modeling errors and error underestimates contribute equally to the increased £ 2 . Therefore 
the x 2 value corresponding to the underestimated errors is 1 + (j 2 - l)/2. Consequently, in order 
to obtain a more realistic error for xl, the formal a x \ is multiplied by this estimated quantity; i.e., 


Table 8. Baseline length scatter ( S B ), improvement (A d B ), and bias (p B ), mm. 


Mapping 

Function 

15—45 

Sr Xl _ 

15-65 

Sb Xl 

15-45 
A 8 b Pb 

15- 

ASb 

-65 

Pb 

1 

49.4 

5.8 

24.6 

8.0 

-40.8 

173 

-20.6 

172 

2 

42.1 

4.2 

24.3 

7.8 

-31.5 

90 

-20.3 

112 

3 

38.0 

3.1 

24.8 

7.9 

-25.8 

26 

-20.9 

9 

4 

27.9 

1.7 

13.4 

2.2 

. . . 

. . . 



5 

27.9 

1.7 

13.4 

2.2 

0.0 

-3 

0.0 

-2 

6 

38.4 

1.7 

19.9 

3.0 

-26.4 

6 

-14.7 

6 

7 

28.1 

1.6 

13.9 

2.3 

-3.3 

20 

-3.7 

15 

8 

29.9 

2.0 

17.3 

3.8 

-10.8 

62 

-10.9 

41 

9 

32.5 

2.3 

17.1 

3.7 

-16.7 

61 

-10.6 

40 

10 

30.3 

2.0 

16.2 

3.3 

-11.8 

58 

-9.1 

45 
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Figure 7a. Normalized of residuals from a linear fit to baseline lengths vs. time for the 
California-Australia DSN baseline. (See text for explanation of error bars.) 


1.3 and 1.6 for the two baselines, respectively. This gives o x i = 0.34 and 0.44, which are the 
error bars used in Figs. 7a and 7b. Since all fits were based on identical data, the statistical 
difference between different solutions is the standard deviation of Xv itself, i. e . , 0.34 or 0.44. 

The differences are due to mapping errors and are seen to be highly significant for the group 
of three Chao functions vs. Lanyi. They are marginally significant for the remaining functions 
on the California-Australia baseline, and > 2(7^ significant for Ifadis, MTT, and NMF vs. Lanyi 
on the Califomia-Spain baseline. The Ifadis, MTT, and NMF mapping functions all worsen the 
scatter by at least 10 mm on both baselines. The CfA-2.2 function resembles Lanyi most closely 
in regard to scatter, and the two functions appear to be nearly identical. The Lanyi map with 
parameters estimated achieves its residual improvement at the expense of increases in baseline 
scatter of 26 and 15 mm. The purpose of this fit is to model any remaining variation in the 
mapping function; indeed, the RMS observable residuals decrease. This improvement appears to 
be achieved by propagating the delay residual errors into systematic baseline errors. The 
mechanism is poorly understood and is currently under further study. Most mapping functions 
also yield substantial baseline length biases relative to Lanyi [1984], ranging from 20 to 61 mm 
for the California-Australia 10,600-km baseline (the Chao tables do very well in this regard, 
only shifting the scale by 26 mm). 
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Figure 7b. Normalized X 1 of residuals from a linear fit to baseline lengths vs. time for the 
California-Spain DSN baseline. (See text for explanation of error bars.) 

5.4 Remarks 

Based on the comparisons discussed in the previous subsections, a number of conclusions 
can be drawn. In the absence of more generalized testing, these should be regarded as limited to 
mapping the tropospheric delays at the three DSCCs. The statistical significance of differences 
among the mapping functions is tied to the residual scatter of these particular DSN VLBI fits. 
First, all three variants of the Chao mappings were very inadequate. They produced post-fit 
residual and baseline scatter values that were inferior by many cm relative to the more modem 
functions. Second, the 1984 Lanyi mapping function is still, by a small margin, equivalent or 
superior to all the newer algorithms developed during the intervening decade. Third, those 
mapping functions that employ limited or no surface meteorological data are either equivalent to 
or slightly worse than the Lanyi function. Some of their deficiencies may be due to the fact that 
their functional forms are based on atmospheric profiles measured at predominantly North 
American sites. Finally, of the functions using minimal or no surface data, either the Ifadis or 
NMF functions produce tropospheric calibration models of DSN VLBI measurements that are 
nearly equal in quality to those given by the Lanyi mapping function on the Califomia-Australia 
baseline, but worse on the Califomia-Spain baseline. 
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SECTION 6 


SELECTION CRITERIA 


There are number of selection criteria that can be used to determine whether or not a new or 
enhanced model should be incorporated into an institutional software set. For choice of a 
candidate tropospheric mapping function, the most immediately obvious selection criteria are 
(not necessarily in order of importance): 1) accuracy; 2) ease of implementation; 3) 
computational complexity, i.e., processing time and memory requirements; 4) model complexity; 
5) ease of use, i.e., are surface weather data necessary or can standard inputs be overridden by 
the user; 6) experience base, or “maturity,” of a particular model; and 7) “tuning” capability, i.e., 
are partials available for estimated and “consider” parameters which can be adjusted. 3 In Table 
9, an attempt is made to grossly categorize each of the candidate mapping functions against these 
seven selection criteria, each of which is discussed in greater detail below. 


Table 9: Cursory rating of candidate tropospheric mapping functions against various selection 
criteria (on a scale of 1 to 4, where 1 is a low rating, and 4 is a high rating). 




Selection Criteria 


Mapping 

Function 

Accuracy 

Ease of 

Implementation 3 

Computational 

Complexity 

Model 

Complexity 

CfA-2.2 

4 (> 5°) 

3 

3 

4 

Lanyi 

4 (> 6°) 

2 

2 

1 

Ifadis 

4 (> 2°) 

3 

3 

4 

MTT 

4 (> 2°) 

3 

3 

4 

NMF 

4(>3”) 

4 

3 

3 


Ease of 
Use 

Maturity 

Tuning 

Capability 


CfA-2.2 

3 

4 

2 


Lanyi 

2 

4 

4 


Ifadis 

3 

4 

2 


MTT 

3 

3 

2 


NMF 

4 

1 

2 



a ODP specific 


3 Consider parameters are an integral part of a “consider state analysis.” a special case of reduced-order filtering theory used in 
orbit determination analysis. Various systematic error sources are treated as unmodeled parameters which are not estimated, but 
whose effects are accounted for (i.e., “considered”) in computing the error covariance of the estimated parameters. 
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6.1 Accuracy 


Although difficult to quantify, each candidate mapping function possesses a range of 
elevation angles and atmospheric profiles that demonstrates superior accuracy characteristics 
over its counterparts. Based on the accuracy assessment presented in the previous section, all of 
the candidate mapping functions demonstrate superior performance characteristics over the 
currently implemented Chao standard mapping tables and approximation formulas. A number of 
studies have been conducted within the past year or two that attempt to quantify the accuracy 
characteristics of various modem tropospheric mapping functions, including those presented in 
this report. For other recent VLBI-based results, see MacMillan and Ma [1994], and for recent 
GPS-based results, see Mendes and Langley [1993]. It should be noted that neither of these 
studies included the NMF mapping functions, as the studies preceded its release. However, the 
recent Mendes-Langley ray trace study does include the NMF functions [Mendes and Langley, 
1994], As this continues to be an active area research, and there is little doubt that by the time 
this report is printed, there will be a number of newly published reports that will give greater 
insight into the performance of today’s highly accurate tropospheric mapping functions. 

6.2 Ease of Implementation 

Computer program mechanizations exist for all of the candidate mapping functions 
presented. Therefore, actual coding of the algorithms is not an implementation concern. The 
issue of accommodating a particular mapping function into the ODP’s infrastructure is, however, 
a concern. Recall that the standard ODP tropospheric calibration model treats the dry and the 
wet components of the tropospheric adjustment separately. Therefore, with the exception of the 
Lanyi mapping function, any of the candidate dry and wet mapping functions could easily be 
incorporated into the ODP tropospheric adjusts model simply by replacing ELVFAC in Eq. (9) 
with the appropriate procedural call. Because all of the candidate mapping functions (with the 
exception of NMF) have provision for specifying real-time surface meteorological data, further 
ODP modifications would have to be made in order to accommodate these additional inputs. 
This would most likely be accomplished in the form of additional inputs on the General INputs 
(GIN) file. Such a requirement is necessary not only for the computations of the actual 
corrections to the modeled observables, but also for the computations of the parameter partials 
that the user may wish to adjust. Therefore, both the EDIT and REGRES program links would 
be impacted. In addition, the Lanyi, MTT, and NMF mapping functions all utilize the site- 
dependent parameters of geodetic station latitude and height. Therefore, procedural calls to 
special utility libraries would have to be made in both program links in order to convert the 
station coordinates given in the conventional spherical or cylindrical coordinates relative to the 
prime meridian, equator, and mean pole of 1903.0 to geodetic station coordinates. Finally, the 
NMF dry mapping function requires the day-of-year (DOY) as a functional input; thus, an 
appropriate time conversion would be necessary to convert the time past epoch from the current 
ODP time format to DOY. 
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6.3 Computational Complexity 

Because high-speed workstation computers are readily accessible in the work center today, 
computational limitations with respect to the different candidate mapping functions are not a 
concern, at least not for tracking and navigation operations at JPL. As an example, the 
processing of several thousand data points that were analyzed in the DSN Mark III VLBI 
accuracy assessment was performed with no appreciable difference between mapping functions 
in terms of processing time or memory requirements. Therefore, no attempt was made to 
benchmark the various candidate mapping functions for time and space complexity on different 
computer platforms. A cursory examination of the processing time requirements for these and 
several other mapping functions (including Chao’s original approximation formula) is presented 
in Mendes and Langley [1993], Results from that study suggested that for certain computational 
platforms, the Lanyi mapping function’s computational requirements were significantly greater 
than the other mapping functions evaluated. 

6.4 Model Complexity 

The level of mathematical sophistication or complexity of a proposed model may be targeted 
as one of the critical criteria of choice over another candidate model. For the majority of analysts 
who would invoke a selectable tropospheric calibration model in the ODP, the level of 
complexity of the underlying model is generally not an issue. Users are typically more 
concerned with the model’s ease of use rather than its mathematical sophistication, as long as the 
model has been proven to yield improved results over the former model. In contrast, if an 
analyst wishes to make modifications or enhancements to a newly incorporated model, the issue 
of model complexity becomes much more important. Upon review of the mathematical 
formulations of the candidate tropospheric mapping functions presented in this report (cf., 
Section 4), it is clear that the Lanyi model is by far the most complex, as it was derived 
analytically. Although the input parameterizations vary significantly between the CfA-2.2, 
Ifadis, and MTT semi-empirical models, all are relatively straightforward and similar in format. 
The NMF model, although somewhat more sophisticated in its parameterizations of the mapping 
coefficients than the other empirically based models, is still relatively straightforward in its 
formulation. 

6.5 Ease of Use 

Another important factor in the selection process of a proposed model upgrade, and 
somewhat related to model complexity, is how easy it is to use. In the current ODP 
implementation of the tropospheric mapping phase, minimal input is required of the user in order 
to utilize the Chao dry (TABDRY) and wet (TAB WET) standard mapping tables for data editing, 
and the Chao dry and wet approximation formulas for the parameter partials. Although the 
capability for overriding the TABDRY and TABWET arrays has existed for over 2 decades, this 
option has yet to be exercised. Because the majority of the candidate mapping functions support 
surface meteorological parameters, there is a trade-off between the potential for greater accuracy 
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and the additional requirement on the user to maintain much greater attention to detail. Standard 
or default values for the surface meteorological parameters could be used for routine 
applications, resulting in less accuracy, but with reduced burden on the user. Because the NMF 
dry and wet mapping functions do not rely on input surface meteorological parameters, they are 
clearly the easiest mapping functions to use, with all inputs transparent to the user. 

6.6 Maturity 

The experience base, or “maturity,” of a particular model must also be considered in the 
evaluation process. The Lanyi mapping function has been used extensively at JPL for precision 
VLBI- and GPS-based analysis since its pre-release in 1983, while the CfA-2.2 dry mapping 
function of Davis et al. has also seen extensive use, principally by the U.S. East Coast and 
international geodetic science activities, since the mid-1980s. The MTT mapping functions, 
apparently intended to supersede the CfA-2.2, were published in 1992 and have also seen 
extensive use in high-precision geodetic applications. The global Ifadis mapping functions, 
although introduced as early as 1986, do not appear to have been as actively “marketed” as some 
of the other candidate mapping functions presented [ McCarthy , 1992]; nevertheless, their 
performance in the precision applications has thus far has been commendable, yielding results 
comparable to or better than the MTT mapping functions. The NMF mapping functions are 
clearly the “new kid on the block,” and as such, do not possess the level of maturity seen for the 
other candidate mapping functions. However, preliminary analysis against ray traces based on 
radiosonde profile data [Mendes and Langley, 1994] and experimental evaluation using VLBI 
measurements (c/, Section 5), suggests accuracy levels as good or better than the Ifadis and 
MTT mapping functions. 

6.7 Tuning Capability 

The final selection criteria suggested herein regards a candidate model’s “tuning” capability, 
i.e., the capability of a particular model to support parameters which can be adjusted in the data 
filtering or estimation process. More specifically, partials are available for parameters which can 
either be estimated or treated as unmodeled “consider” parameters. The most commonly utilized 
error model parameter partials for the tropospheric calibration model are with respect to zenith 
dry and wet delays [cf, Eqs. (21) and (22)]. In an effort to further improve the parameter 
estimation process, the mapping coefficients given by the CfA-2.2, Ifadis, MTT, and NMF 
mapping functions could be promoted to estimable parameters, as was done for Chao’s original 
semi-empirical dry mapping function in JPL’s VLBI parameter estimation software set, 
MODEST [ Sovers and Jacobs, 1994], Currently, only the Lanyi mapping function supports an 
additional tuning parameter to account for the absence of real-time surface meteorological data. 
Recall that in the reduction of DSN Mark III VLBI measurements (discussed in Section 5), use 
of this additional tuning parameter resulted in improved post-fit residuals, albeit at the expense of 
increased baseline scatter. 
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SECTION 7 


SUMMARY AND CONCLUSIONS 


The current implementation of the ODP’s standard tropospheric calibration model was 
reviewed in extensive detail. The standard model is an adaptation of the seasonal model 
developed by C. C. Chao in the early 1970s to support the Mariner Mars 1971 program, and has 
seen only modest changes since it was incorporated into the ODP. As part of a model upgrade 
effort, this report reviewed the format and content of the zenith delay tropospheric calibrations 
supplied by the DSN’s TSAC activity, as well as the tropospheric mapping process, whereby the 
zenith delay calibrations are projected along the station-to-spacecraft line of sight. All 
mathematical descriptions of the standard model were presented, along with ODP program link 
interface requirements. Proposed improvements to the current model were described for both the 
zenith delay calibration and zenith delay mapping processes. 

The proposed improvements for the zenith delay calibration process are currently under 
review by TSAC, and include calibrations based on real-time surface meteorological 
measurements, as well as calibrations derived from GPS ground-based observations. The 
principal focus of this report was to review a set of relatively new models in an effort to 
modernize the current zenith delay mapping process. These new models, or “tropospheric 
mapping functions,” are accurate over a wide range of observer-to-radio source elevation angles. 
Moreover, nearly all of the newer mapping functions have been used successfully in a variety of 
high-precision geodetic experiments. Candidate mapping functions of Lanyi, Davis et al. (CfA- 
2.2), Ifadis (global solution model), Herring (MTT), and the recent Niell (NMF) mapping 
functions were described. All of the candidate mapping functions demonstrated superior 
accuracy characteristics over the currently implemented Chao model, which takes the form of 
standard mapping tables derived from early ray trace calculations and simple approximation 
formulas to perform the mapping. These performance results were based on a statistical analysis 
of DSN Mark III VLBI measurements acquired over a 5-year period from 1988 to 1993. 

Clearly, no one “best” tropospheric mapping function exists for every application and all 
ranges of elevation angles; however, based on the comparative survey presented, the authors 
recommend that the Lanyi and NMF mapping functions be incorporated into the ODP as user- 
selectable options, in addition to the currently implemented Chao model. 

Justification for this recommendation follows: 

1) Although the Lanyi model would admittedly be the most difficult to incorporate into the 
ODP’s current media calibration infrastructure, it has the advantage of being adaptable to a wide 
variety of atmospheric conditions, given the provisions to input real-time weather data. In the 
accuracy analysis presented, the Lanyi mapping function was seen to be superior over its 
counterparts-the CfA-2.2, Ifadis, MTT, and NMF functions-at 6“ elevation or higher. 
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Furthermore, the CfA-2.2 lacks a companion model for the wet atmospheric component. A 
cursory examination was made of the sensitivity of the Lanyi mapping function to variations in 
surface meteorological data. However, the level of accuracy given by the DSN’s surface 
instruments has yet to be determined and will require further investigation; standard (average) 
values for the weather data can be used in the interim that will still yield high-fidelity results. 
Another argument for implementing the Lanyi model is that line-of-sight tropospheric delay 
predictions as measured by the GPS ground system are typically mapped to the zenith direction 
in time series format using Lanyi’s mapping function. Therefore, the Lanyi model would have to 
be an available option in the ODP in order to support any future GPS-based zenith delay 
calibrations provided by TSAC; otherwise, post-fit results will be biased if a different mapping 
function is used to “de-map” the zenith delays. The final argument given here for implementing 
the Lanyi mapping function is that for most applications supported by the DSN, the Lanyi 
mapping function is expected to perform as well or better than the other candidate models at 6° 
elevation angle and above; recall that operational constraints of the DSN radio antennas usually 
preclude tracking below 6°. For non-DSN applications that require tracking significantly below 
this level, the Ifadis, MTT, or NMF mapping functions are recommended. 

2) The NMF model, although lacking maturity, has the advantage of perfonning as well as or 
better than its counterparts, such as the Ifadis or MTT models, without requiring surface 
meteorological input data. Therefore, ODP users would not be burdened with specific model 
interfaces. It appears that the NMF mapping functions will gain popularity over time, since they 
are easy to use and have been shown to be highly accurate over a wide range of elevation angles 
in certain experimental applications. In addition, the NMF mapping functions utilize the same 
form of the Marini continued fraction as Ifadis and MTT (truncated at three terms in the 
expansion), the principal difference being in the parameterizations of the mapping function 
coefficients. Thus, an argument can be made that the NMF model will eventually supersede the 
Ifadis and MTT models. 

3) The Chao model should remain as the default option, at least in the near term. It would be 
imprudent to suggest replacing the Chao model in its entirety, as orbit determination analysts 
frequently attempt to reconstruct “old” data and validate anomalies encountered during past 
mission operations. In addition, the ODP is a “flight-certified” software set; thus, any proposed 
modeling improvements would have to be rigorously tested against the original standard model 
before they could be accepted as a new standard. 

In closing, recall that all of the candidate mapping functions were derived under the 
assumption of a spherically symmetric atmosphere. Except for ultra-high-precision applications, 
this azimuthal asymmetry assumption is typically adequate. Consequently, for routine support of 
flight operations, such a level of model sophistication does not yet seem warranted, at least not 
for ODP implementation. The interested reader is referred to publications by Gardner [1977] 
and Wilcox [1994] for examples of novel solutions to the azimuthal asymmetry problem. 
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APPENDIX 
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